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Abstract 

The autocorrelation function C^^e(A(^, Ae) = {6g{ip, e) 5g{ip + A(^, e + Ae)) 
{if and £ are rescaled magnetic flux and energy) for the magnetoconductance 
of a balhstic chaotic quantum dot is calculated in the framework of the su- 
persymmetric non-linear cj-model. The Hamiltonian of the quantum dot is 
modelled by a Gaussian random matrix. The particular form of the symme- 
try breaking matrix is found to be relevant for the autocorrelation function 
but not for the average conductance. Our results are valid for the complete 
crossover from orthogonal to unitary symmetry and their relation with semi- 
classical theory and an 5-matrix Brownian motion ensemble is discussed. 
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Typeset using REVTgX 



GaAs/Ala;Gai_2^As heterostructures provide useful experimental realizations of two- 
dimensional ballistic cavities known as quantum dots. Measuring the conductance of a quan- 
tum dot connected to electron reservoirs, one can study the quantum behavior of classically 
chaotic billiards. Apart from the channel number, the magnetic field B and the gate volt- 
age E are the two important adjustable parameters. In experiments, the sample dependent 
conductance is measured as a function of these parameters. Theoretical approaches 
consider the weak localization peak (by calculating the average conductance as a function 
of magnetic field) and the conductance autocorrelation with respect to magnetic field and 
gate voltage. Applying an ergodicity argument, one may compare these quantities with 
experiment. 

Chaotic ballistic systems were treated numerically and semiclassically in 0-^] , the latter 
predicting an algebraic decay of the correlation functions and the weak localization peak. 
A random matrix approach for the unitary scattering matrix in terms of the circular Brow- 
nian motion ensemble yields an exponential decay of these quantities as a function 
of a fictitious Brownian motion time t. It is not obvious how to relate these results with 
the semiclassical theory. Furthermore, it is not entirely clear if the problematic diagonal 
approximation in the Gutzwiller trace formula, used in Refs. [§,0, plays an important role 
here. 

Recently, the supersymmetric non-hnear a-model has been applied to quantum dots 
to calculate the weak localization peak, describing the crossover from orthogonal to unitary 
symmetry for the average conductance. In the unitary symmetry class where the magnetic 
field is sufficiently strong to completely break the time reversal invariance, the conductance 



autocorrelation function has been calculated [10] in the limit of many channels. 



Corresponding perturbative treatments of the a-model which are predominantly con- 
cerned with the metallic diffusive case can be found in Ref. |TT[] for the magnetic field 



correlations and in Ref. [T^ for the energy correlations, respectively. 

The purpose of this letter is to generalize the treatment of Ref. to the complete 
crossover from orthogonal to unitary symmetry. In addition, we point out that the precise 
form of the symmetry breaking matrix is relevant for the autocorrelation function but not for 
the weak localization peak. This observation is very important to understand the Brownian 
motion results 

We are interested in the generic chaotic features of the quantum dot and model the 
Hamiltonian as a Gaussian random matrix of the form 

^a=i^SE + «('^^SE + ^^) (1) 

where Hqq-^ are (independent) real symmetric iVx A^-random matrices given by the Gaussian 



orthogonal ensemble |jT3| with the variance X'^/N of the nondiagonal elements and A is a 
real antisymmetric random matrix whose independent elements are normal distributed with 
zero mean and variance X'^/N. The parameter a (which is proportional to the magnetic fiux 
penetrating the quantum dot) describes the strength of the symmetry breaking. Concerning 
the parameter k, which distinguishes between different types of symmetry breaking, we are 
mostly interested in the particular cases of a purely imaginary antisymmetric symmetry 
breaking (k = 0) and a hermitian symmetry breaking {k = 1). The first case yields the well 



known Pandey-Mehta Hamiltonian |14] which is believed to correspond to the application 
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of a small magnetic field [0,0]. The other case is rather directly related [0 to the circular 
Brownian motion ensemble. 

In the large limit, A is expressed in terms of the average level spacing Aq (at energy E = 
0) |13| by A = N Aq/tt and the proper crossover-parameter [|l^] to describe the symmetry 



breaking is just v Na. Keeping this quantity finite as ^ oo the value of k does not affect 
S-matrix averages because the additional contribution of Hqq^ can be taken into account by 
an infinitesimal rescaling of the variance of -f^QOE- ^^e other hand, correlations between 
different values of a depend on the choice of k. 

We follow the treatment of Refs. and use for the M x M-scattering matrix the 

expression S = 1 — 2'KiW\E — + mWW'^)~^W where 1^ is a A^ x M-matrix describing 
the coupling of the scattering channels with the states of the quantum dot. As in Ref. 0, 
we use the choice of W that corresponds to the ideal coupling of the wires characterized 
PJT6[| by a vanishing average (S*) = and leading to equivalent scattering channels. Then, 



the distribution of the scattering matrix (for a = 0) is |1T7|JT8[| equivalent to the circular 
orthogonal ensemble |19[. We assume the dimension M of the 5- matrix to be even and 



express the conductance by the Landauer formula g = Tr(S'|2'S';^2) where S12 is the M/2- 
dimensional (1, 2)-block of S. 

We apply the supersymmetric approach as described in Refs. [p!6| , p0| -p^ to calculate the 
conductance autocorrelation function at two different energies E ± AE/2 and two different 
values a ± Aa/2 of the symmetry breaking parameter. We assume that in the large A^- 
limit the quantities A^a^, A^Aa^, E/Aq, and AE/Aq remain finite. The product of the 
two conductance contributions is expressed as a Gaussian integral over a 16A^-dimensional 
supervector with the 16A^-dimensional supermatrix 

H = A^E + ^j:,~[HSl)E + i(^ + ^^3)i'^Hg!,^ + tAT3))+t'KAWW^^ . (2) 



We adopt here the notational conventions of Ref. ||2T|. The supermatrices A, and S3 are 
diagonal matrices with an equal number of diagonal entries +1 and —1, respectively, defining 
different types of gradings. A describes the grading imposed by the advanced and retarded 
Green's functions and the rs-grading is determined by the time reversal transformation. S3 
corresponds to the additional grading for the two conductance contributions at different 
magnetic fields and energies. As in we denote by Lg the supermatrix describing the 



decomposition in bosonic and fermionic sub-blocks. 

Performing the usual steps p!6| , |20| -p2| and using the rescaled quantities (p = {SNa"^ /{M + 
1))^/^ and e = 27iE/[(M + l)Ao] for the symmetry breaking parameter and the energy, 
respectively, we arrive at the Q-integral 

g{if - ^Acp, e - \Ae) g{ip + ^Acp, e + \Ae)) (3) 



2\ 2 



dQ n [Str((l + gA)-iQAJ, 



u=l,2 

with the action 

C{Q) = ^Str ln(l + QA) + (-4z A^ Str (gS3) + ^ Str ((Q7.)')) (4) 



1^=1,2 
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and the matrices 71 = + ^AipEs^T^, 72 = ^AipE^, Ji^2 = (1 + Lg){l ± E3)(Ti(A)/4 
((Ti(A) is the first Pauh matrix in the A-grading). The integration variable Q is a 16 x 16- 
supermatrix which belongs to the coset-space defining the non-linear cr-model [|l|,2C,2l|] for 
the orthogonal symmetry class. 

Now, we consider the limit M ^ 1 and perform an expansion in terms of the small 
parameter x = {1 + M)^^. In typical experiments with ballistic quantum dots, the number 
of channels is about 3-6 and correspondingly x = 1/7-1/13. Using the well known square root 



parametrization [|TB],p2|, the Q-matrix can be written as Q = T ^AT, T = y/l + D^ + D. D 



is a supermatrix which has only non diagonal entries in the A-grading: Dn = D22 = and 
D21 = LgDi2. The integration measure [^] is dQ = exp ^Str ln(l -|- -D^)) dj2{Di2) where 
dn{Di2) is the flat measure of the 8x8 supermatrix D12. The Jacobian factor of this measure 
modifies the action in Eq. (^, i. e. in the first contribution M is be replaced by M -|- 1. 
The Q-integral can then be evaluated perturbatively p2[ where the (modified) action has to 
be expanded in even powers of D, C{D) = C-2{D) + C4,{D) + Cq{D) + ■ ■ ■. The expansion of 
the source term contributions and of the exponential of Ca,{D) + Cq{D) up to three leading 
orders in x yields a Gaussian integral with the action €2(0). The supermatrix D12 has to 
be decomposed into the 16 blocks corresponding to the T3- and Sa-grading. Each of these 
blocks is associated with a particular propagator determined by £2(0). Due to the large 
number of different terms in the higher order supertraces, we applied a computer program, 
which was developed for this purpose in C-I--I-, for the final evaluation of the integral. To 
get the autocorrelation, one has to subtract from (^) the product of the two independently 
averaged conductance contributions which have been calculated similarly. 

Before we discuss the final result, we mention that the action (^) also describes diffusive 
metallic systems [pO| , ^ provided that all relevant energy scales, in particular the level 



broadening M ■ Aq, are smaller than the Thouless energy. The perturbation theory applied 
here corresponds exactly pO| to the standard diagrammatic approach in terms of diffuson 



and cooperon modes. The relevant contributions for the conductance autocorrelation are 



the diagrams consisting of two connected conductance loops. In the cr-model, they 
correspond |2^JTl|Jl^ to the contributions of D12 that are non diagonal in the Sa-grading 
and are associated with two propagators: the cooperon mode Pc = [(M-|-l)(l-|-v9^ + ^Av9^-|- 
iAe)]~^ (non diagonal in the rs-grading) and the diffuson mode Pd = [{M+l){l + ^^^Aip'^ + 
iAe)]~^ (diagonal in the ra-grading). The conductance autocorrelation is estimated as 
~ M^(|PdP -|- I -Pep) which is indeed confirmed by the explicit computer calculation. We 
obtain (after the shifts ip ip+^Aip, e e+^Ae) for the average conductance {g{f, e)) and 
the autocorrelation function C^^e(A(y9, Ae) = {Sg{ip, e) Sg{ip + Aip, e + Ae)) {6g = g — {g)) 
the results 

^^^^'^)^ = T-4(MTT)rT^ ' 

C^,,(A(^, Ae) = ^ (- ^—^ -2 + 2 ) > (6) 

16V(l + (<^ + iA(^)2 + ^jA(^2) +^^2 ^^^2 J 



which are correct up to terms of order (9(x^) for (|^) and of order 0{x) for (|^). 

The quantities v? = (8iVaV(M-M))^/^ and Ae = 27iE /[{M+l)Ao\, given in terms of the 
original model parameters, are of course related to the corresponding physical quantities in 
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a ballistic quantum dot. The energy correlations decay on a scale 7 = Ao(M + l)/7r which is 
just the level broadening (or inverse life time) due to the coupling with the channels . Aq 
has to be adjusted to the level spacing in the quantum dot at the Fermi energy. The relation 
between the parameter \fNa and the magnetic flux $ applied on the dot was studied in 



||15|,[7[| and is given by a = const. \JhvF / {L^o) ($/$o) where vp is the Fermi velocity, 
L a typical diameter of the quantum dot and $0 = h/e the flux quantum. The numerical 
constant is of the order of unity and depends on geometrical details. For a particular model, a 



circle (radius L) with a very rough surface, it takes the value y4/3 ~ 1.15 0. To estabhsh 
the connection with experiment one should therefore set = const. (Shvp/L'jy^'^ (^/^q) 
and Ae = AE/{2'y) in Eqs. (^ and (^. The scale dependence on the channel number M 
(via 7) can be understood in terms of semiclassical trajectories. With increasing number of 
channels (decreasing life time h/j) the typical length of the trajectories (and therefore the 
enclosed area) decreases so that the critical value for the energy (or the magnetic flux) to 
break the phase-coherence of the electrons is increased. 

In the limiting cases of orthogonal (/3 = 1, 9? = 0) and unitary (/5 = 2, — > 00) symmetry, 
the correlation function becomes a squared Lorentzian in Aip and a simple Lorentzian in Ae 

C,,,(A^,A5) = i-- / , . (7) 

8P (i + ^Aip^) + A£2 

For K = 0, (3 = 2 this expression is consistent with the semiclassical approach and the 
result of Ref . [l^ , provided the scales of the flux are properly translated. The prefactor in 
(0) exhibits the /5-dependence of the universal conductance fluctuations for a ballistic cavity 
obtained by the random matrix approach for the scattering matrix |^,25|. The expression 



(^) for the average conductance is included here for completeness and agrees very well with 
the (numerical flt of the) precise result of Ref. 

The semiclassical calculations of Refs. [^|| do not reproduce the correct amplitude but 
they predict the same functional dependence of the weak localization peak and the autocor- 
relation function on the magnetic fleld as (^ and (0), provided the imaginary antisymmetric 
symmetry breaking (k = 0) is considered. We also recover the simple Lorentzian for the 
energy correlation Q. 



In Refs. p^ , p5| the ^-matrix of the quantum dot was described by one of the circular 
ensembles for orthogonal or unitary symmetry. The crossover between these cases 
can be modelled by a so-called Brownian motion ensemble where the S'-distribution 
depends on a flctitious time t corresponding for t = to the circular orthogonal ensemble 
and "diffusing" for t — > 00 to the circular unitary ensemble as stationary distribution. The 
average conductance and the autocorrelation function obtained in this approach are given 
by 01: 

(git)) = ^ - e-*/*^ , {6g{t) 6g{t + At)) = ^ e"^*/*^ (l + e'^*/*^) + 0{M'') 

(8) 

where tc is the critical time that determines the crossover scale. The comparison of the 
average conductance expressions (H) and (13) suggests the identiflcation t = tcln(l + v^^) 
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between the Brownian motion time t and the parameter ip. Then, the conductance fluctua- 
tions given by (Ay9 = Ae = 0) and (|^) {At = 0) coincide. Concerning the correlations, 
we obtain from (|^) (for the limiting cases if = or (p ^ oo) in lowest order in Ay9^: 
C^^£(Ay9, 0) ~ g^(l — |(1 + K^)Ay9^), which coincides with the Brownian motion expression 
(with At/tc — A(y9^) only for the case of hermitian symmetry breaking (k = 1). It is indeed 
possible to relate the Brownian motion approach directly (for small ip, Aip, large M) 
to the model considered here (k = 1), giving the correct identification Aip"^ = At/tc- In 
the crossover regime ip ^ 1, this identification and the mapping onto the Brownian motion 
approach do not work since the correlation function (^ is not even in Aip. 

In summary, we have calculated the conductance autocorrelation function (with respect 
to magnetic field and energy) of a chaotic quantum dot for the complete crossover from 
orthogonal to unitary symmetry. The perturbative approach seems to work rather well, 
since our expression (|^) for the weak localization peak agrees very well with the precise 
result of Ref. In addition, we have found that different types of symmetry breaking lead 
to different scales in the autocorrelation function explaining the inconsistencies between 
semiclassical theory (/t = 0) and the Brownian motion approach (k = 1). 

* * * 

The author thanks J.-L. Pichard, A. Miiller-Groeling, H. A. Weidenmiiller, and C. W. J. 
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